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Abstract 

We discuss the nuclear dependence of •i/' and ^' production in hadron-nucleus 
interactions as a function of longitudinal momentum fraction xp- Nuclear ef- 
fects such as final-state absorption, interactions with comovers, shadowing of 
parton distributions, energy loss, and intrinsic heavy-quark components are 
described separately and incorporated into the model which is then compared 
to the preliminary E866 data. The resulting nuclear dependence of Drell-Yan 
production at 800 GeV and proposed measurements of V', V'' and Drell-Yan 
production at 120 GeV are also calculated. 
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1 Introduction 



The factorization theorem of perturbative QCD [0], separates the perturbatively- 
calculable short- distance quark and gluon interactions from the nonperturbative dy- 
namics underlying the parton distribution functions in the hadron. The effectiveness 
of the factorization theorem in nuclear targets with mass number A can be obtained 
by a comparison of the perturbative production cross sections of hard processes in 
nuclei to those in a free proton, or, since nuclear isospin effects are generally small, 
a nucleon. The dependence of particle production on atomic mass number A is con- 
ventionally parameterized by a power law as 0, H, ||, ||, |^, ||, ^ 

(^pA = CTpNA" , (1) 

where apA and apN are the integrated particle production cross sections in proton- 
nucleus and proton-nucleon interactions respectively. If factorization is satisfied, then 
particle production should be independent of the presence of nuclear matter and 
apA would grow linearly with A, implying a = 1. Drell-Yan production, integrated 
over all kinematic variables, shows this linear growth to rather high precision 0. A 
number of experiments have measured a less than linear A dependence for ip and 
ip' production ||, ||, ||, ^ |], ^. Typical values of the exponent a in Eq. (|l|) are 
between 0.9 and 1. This nonlinear growth of the ip cross section with A has been 
used to determine an effective nuclear absorption cross section 0, |10|, |lT| under the 



assumption that the deviation of a from unity is solely due to ip dissociation by 
nucleons. However, attributing the entirety of the integrated nuclear dependence to 
final-state absorption neglects other possible contributions, perhaps resulting in an 
overestimate of the nuclear absorption cross section. This paper identifies a number 
of possible nuclear effects on ip and, consequently, Drell-Yan production and examines 
their contributions to Eq. (0). 

Any dependence on the kinematic variables such as projectile energy or longitu- 
dinal momentum fraction, xp, would reveal the importance of going beyond a simple 
A scaling for production and a constant absorption cross section for ip production. 
Indeed, it has long been known that in quarkonium production a decreases as a func- 
tion of P, H, I^. There are a number of effects which could contribute to the xp 
dependence. The nuclear parton densities are systematically different from those in 
a deuteron or a free proton as a function of parton momentum fraction x |jl2[. Such 



alterations of the parton densities are universal because they are independent of the 
final state, affecting both ip and Drell-Yan production. Another possible universal 
component is energy loss by the incoming parton as it traverses the nucleus [ll3| , |14|. 



The remaining mechanisms primarily affect ip and ip' production. Absorption of the 
produced ip or cc state by interactions with nucleons and/or produced particles has 
been claimed to be responsible for all ip suppression in nuclear collisions, at least up to 
S+U interactions [jlO|, [11], |l^, |l6l. The importance of these absorption effects depends 
on the production mechanism and the magnitude of the interaction cross sections. 
Energy loss by the final-state color octet cc has also been suggested [0. Finally, the 



presence of intrinsic heavy-quark component of the projectile wavefunction are also 
considered [|1^, ^ . 
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We study the xp dependence using a two-component model employing concepts 
developed in Ref. p^. The first component, based on perturbative QCD [^], is a 
hard-scattering approach that would yield an approximately linear A dependence, 
as in dilepton production by the Drell-Yan mechanism. The A dependent effects 
associated with hard scattering are final-state interactions, nuclear shadowing in the 
target and energy loss in the projectile. The second component of the xp dependence 
arises from intrinsic cc pairs in the projectile wavefunction [1^, |T^. Since the charm 
quark mass is large, these intrinsic heavy quark pairs carry a significant fraction of the 
longitudinal momentum and contribute at large xp whereas perturbative production 
decreases strongly with xp- The light spectator quarks in the intrinsic cc state interact 
on the nuclear surface, leading to an approximate A?^^ dependence |jl9| . 

Such a separation of production mechanisms was first proposed by the NA3 col- 
laboration when they divided their data into hard and diffractive components so 
that 

dxp dxp dxp 

The 'hard' component, ah, includes nuclear shadowing and parton energy loss which 
can alter its effective A dependence. Final state dissociation of the ip or cc state by 
nucleons and secondaries, which does not affect the parton densities, also contributes 
to the effective exponent a'. The 'diffractive' component, ad, is assumed to be due 
to intrinsic charm and only contributes significantly at xp > 0.25. The NA3 Col- 
laboration found a' = 0.97 and P = 0.71 for proton projectiles 0]. Taken together, 
these components give the effective a of Eq. (|I|). We will discuss each mechanism in 
detail in the following sections. Since Drell-Yan production would also be affected by 
nuclear shadowing in the target and energy loss by the projectile partons, its effec- 
tive A dependence should also depend on xp. Thus we will include the Drell-Yan A 
dependence in our study. 

We will calculate a{xp) for each effect individually and then compare the shapes of 
a{xp) with the preliminary E866 ip and ip' 800 GeV pA data |0, shown in Fig. |1|, when 
all the mechanisms are combined. For the ip a 0.94 until xp ~ 0.25, decreasing to 
a ~ 0.7 at large xp and for the ip' a ^ 0.91, effectively constant, until xp ~ 0.5. The 
two values of a are essentially compatible with each other within the experimental 
uncertainties. No decrease at negative xp is observed, contrary to previous results 
1^, |[. The drop to a ~ 0.7 at xp ~ 0.9 is similar to previous measurements p, ^. 
Since the E866 targets are tungsten, W {A = 184), and beryllium. Be {A = 9), we 
calculate the cross section per nucleon for each target according to Eq. (^, obtaining 
from Eq. (|l]), 

X _ . , ln[{dapw/Awdxp)/{dapBc/ABedxp)] 

a[xp) — i H — — — — — r . [6) 

\n{A-w/ABe) 

We first discuss and Drell-Yan production in QCD and then each nuclear ef- 
fect in turn. Quarkonium production by color evaporation and in non-relativistic 
QCD is discussed in Section 2 and dilepton production by the Drell-Yan mechanism 
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Figure 1: The preliminary E866 data for the ip and ip' ^ dependence as a function 
of Xp- 

is briefly touched upon in Section 3. Three different models of nuclear absorption 
are discussed in Section 4. The first two, absorption of pure color octet and color 
singlet states respectively, are used in conjunction with the color evaporation model 
of quarkonium production. The last, a combination of octet and singlet absorption, is 
coupled to quarkonium production in non-relativistic QCD. A discussion of quarko- 
nium dissociation by comoving secondaries is presented in Section 5. Three different 
parameterizations of nuclear shadowing are described in Section 6. Several models 
of energy loss are discussed in Section 7. The intrinsic charm model is introduced in 
Section 8. The combined results are given in Section 9, along with predictions for 
projected lower energy measurements at 120 GeV. A nontrivial combination of effects 
is required to understand the ip and ip' data. 

2 Charmonium Production 

There are two basic models of quarkonium hadropro duct ion that have enjoyed con- 
siderable phenomenological success. The first, the color evaporation model, treats 
all charmonium production identically to cc production below the DD threshold. 
The more recent non-relativistic QCD approach involves an expansion of quarkonium 
production in powers of v, the relative Q-Q velocity within the bound state. Each 
approach will be described in turn and the xp distributions in pp interactions will 
be presented to provide a basis for understanding dat/dxp, Eq. (|^), before nuclear 
effects are included. We will also show the relative contributions from gg fusion, qq 
annihilation, and gq scattering. 
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2.1 The Color Evaporation Model 



In the color evaporation model, CEM, quarkonium production is treated identically 
to open heavy quark production except that the invariant mass of the heavy quark 
pair is restricted to be less than twice the mass of the lightest meson that can be 
formed with one heavy constituent quark. For charmonium the upper limit on the 
cc pair mass is then 2m£,. The hadroproduction of heavy quarks at leading order 
(LO) in perturbative QCD is the sum of contributions from qq annihilation and gg 
fusion. The hadroproduction cross section is a convolution of the qq and gg partonic 
cross sections with the parton densities in projectile A and target B. If xp is the 
cc longitudinal momentum fraction in the AB center-of-mass frame and -y/s is the 
center-of-mass energy of a nucleon-nucleon collision, the cross section for production 



of free cc pairs with mass m is |22 



da" 



dxpdm? 



/ dxidx2 6{xiX2S — rn?) 6{xf — Xi + X2)Hab{xi-, X2; m? 
Jo 

H AB (x 01, X 02; m"^) 



(4) 
(5) 



where xi and X2 are the fractions of the hadron momentum carried by the projectile 
and target partons respectively. After integration over the delta functions in Eq. (^), 
3^01,02 = ^(ixi? + \Jx\ + Am? / s). The convolution of the partonic cross sections and 
the parton densities is 



Hab{xi, X2] m^) = fg{xi, w?)ff{x2, m^)agg{m^) 

+ E [/,^(xi,m2)/i^(x2,m2) + 4^(xi,m2)/f(a;2,m2)]Mm2) 



(6) 



q=u,d,s 



where the parton densities fi{x, m?) are evaluated at momentum fraction x and scale 
w? = X1X2S and m is the invariant mass of the cc pair. The sum over q includes only 



light quark flavors. The LO partonic cross sections are p2 



agg{m^) 
agg{m'^) 





(7) 



where A = y 1 — Am1/w?. 

The LO charmonium cross section for charmonium state z, 5"j, is then obtained 
by integrating the free cc cross section over the pair mass from the cc production 
threshold, 2mc, to the open charm threshold, 2mj:) = 3.74 GeV. Then 



ddi 



2F,; 



2m r 



m dm ■ 



da' 



dxp ' J2mc " dxpdm^ ' 
where is the fraction of cr^^ that produces the final-state cc resonance. 



(9) 
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The CEM assumes that the quarkonium dynamics are identical to those of low in- 
variant mass cc pairs. The hadronization of the charmonium states from the cc pairs 
is nonperturbative, involving the emission of one or more soft gluons. A different 
nonperturbative matrix element is needed for the direct production of each charmo- 
nium state. Each nonperturbative matrix element is represented by a single universal 
factor Fi which depends on the charm quark mass, rric, the scale of a^, /i, and the 
parton densities^. In our calculations with the CEM, we use the leading order MRST 



LO parton distributions |]23|. This set, more recent than the GRV 94 LO densities, 
has a low initial Q^, Qq = 1 GeV. Once has been determined for each state, e.g. ip, 
i'' or Xcj, the model successfully predicts the energy and momentum dependencies. 
We note that F^ includes both direct ip production and indirect production through 
radiative decays of the Xcj states and hadronic ip' decays. 

Since Fi must be a constant for the model to have any predictive power, the relative 
differential and integrated quarkonium production rates should be independent of 
projectile, target, and energy. This appears to be true for the charmonium production 
ratios E j Xcj/^ ~ 0.4 and ip'/^ ^ O.U^,^^^^. See Ref. [gg] for more 
details. 

The next-to-leading order (NLO) quarkonium production cross section in the CEM 
|2T[1 was calculated using the QQ production code of Ref. with the mass cut in 



Eq. @. When the NLO contribution is included, the px dependence of ip production 
at the Tevatron has been shown to agree with the CEM calculations [Q. The LO 
and NLO calculations agree equally well with the energy and xp dependent data if 
Fj^^ is defined as Fj^^'^ multiplied by a theoretical K factor, the ratio of the NLO to 



LO cross sections 21 



Figure ^ shows the forward Xp distributions for ip production in pp collision^ in 
the CEM at 800 GeV and 120 GeV using the MRST LO Q parton distributions. 
The ip' distributions are identical except for the relative fraction of ip' production 
below the DD threshold and are thus not shown. Note that at large xp, xp > 0.6 at 
800 GeV and > 0.5 at 120 GeV, qq annihilation is the most important contribution 
to the cross section. 



2.2 Quarkonium Production in Non-Relativist ic QCD 

An alternative model of quarkonium production, the color singlet model [Q, pre- 
dicted that high pp ip production would be dominated by Xcj decays. It also predicted 
that direct ip and ip' production would be rare because a hard gluon emission was 
required to make a color singlet ^5*1 state on a perturbative timescale. On the other 
hand, the CEM is an average over the color and spin of the produced cc pair and 
cannot make such predictions. Soon after the high pp Tevatron data I^S] made clear 



that the hard gluon emission constraint in the color singlet model severely under- 



^The mass and scale parameters are rric — 1.2 GeV and fj, ~ 2toc for the MRST LO |23|, CTEQ 
4L Q, and CTEQ 3L l^j parton distributions and mc = 1.3 GeV and /i = nic with the GRV 
94 LO parton distributions. At next-to-leading order, the fraction F^p is 2.54% for all parton 
distributions ||2l]] . 

■^The xp distributions are symmetric around xp ^ in pp production. 
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Figure 2: The ip xp distributions at (a) 800 GeV and (b) 120 GeV in the GEM. The 
contributions from gg fusion (dashed) and qq annihilation (dot-dashed) are given 
along with the total (solid). 

predicted direct ip and ip' production, the non-relativistic QGD, NRQGD, approach 
to quarkonium production was formulated p6|- This approach does not restrict the 



angular momentum or color of the quarkonium state to only the leading singlet state. 
For example, the final-state ip may produced as a ^Pq color octet state which becomes 
a ijj through nonperturbative soft gluon emissions. Thus the NRQGD model is similar 
in spirit to the GEM albeit with more nonperturbative parameters, as we will see. 
The xp distribution of a charmonium state, C, in NRQGD is 



dxp 
d{i3 - C) 



: J2 dxidx26{xF - xi + X2) fi (xi, iJ,"^)/ f{x2, fi'^)d-iij C) (10) 



where the C production cross section, a{ij C), is the product of expansion co- 
efficients, C^QQ^j^y calculated perturbatively in powers of Q;s(/i^) and nonperturbative 

parameters, {O^), describing the hadronization of the charmonium state. We use the 
parameters determined by Beneke and Rothstein for fixed-target hadropro duct ion of 
charmonium with rric = 1.5 GeV and n = 2mc and the GTEQ 3L parton densities 
The total ip xp distribution includes radiative decays of the XcJ states and 



37 



hadronic decays of the ip', 



da^ 
dxp 



d<yf 

dxp 



dxp 



da. 



dxf 



:i2) 



In contrast, in the GEM, the xp distributions of all states are assumed to be the 
same. Thus in Eq. (^ implicitly includes the Xcj and ip' decay contributions 
given exphcitly in Eq. (p!2D . 
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For completeness, we now present the individual contributions to ip' and ip pro- 
duction also given in Ref. To simplify the cross sections, we define the coefficients 
proportional to and as 



5, 



(2mc 



(13) 



and denote the delta and theta functions of argument X1X2 — 'imi/s by 5^ and O^- 



Direct production has only contributions from gg fusion and qq annihilation |^ 
as in the CEM, 



^{99 i^) 



e; on 



I - z"^ + 2z\nz 



a{qq V) 



16 
27 



1 — z"^ — 2zln z 
B26,{OtCS,)) , 



(14) 
(15) 



where z = 4m^/(xiX2s). Note that in the terms proportional to 6^, the momentum 
fractions xi and X2 are the same as xio and X20 in the CEM for m = 2mc. The ip' cross 
sections are identical except for the values of the matrix elements Ag, (Oi(^S'i)), and 
(08(^5*1)). The nonperturbative matrix elements for ip and ip' production in Eqs. (|14D 
and (P) are 

A8(^') = 0.0052 GeV^ , 
(Of (3^1)) = 0.76 GeV^ , 
{Oti^Si)) = 0.0046 GeV^ . 



Asi^p) = 0.03 GeV^ , 
(Ofi^Si)) = 1.16 GeV^ , 
(Oti^Si)) = 0.0066 GeV^ 



(16) 

At 800 GeV ~ 63% of the ip' cross section is color octet while 85% of direct ip 
production comes from the octet contribution. This is because As{ip) ~ 5.8As{ip'), 
see Eq. (|TBp. 

Most of the octet contribution to ip production is from direct production. The 
singlet contribution, ^ 40% of the total production cross section at 800 GeV, is due 
to the Xci- The Xcj cross sections from Ref. |3^ are given below according to their 
production mechanisms. The following symmetry relations. 



(2J+l)(Of"(3po)) 



(17) 
(18) 



are particularly useful for calculating the Xcj contributions to the total ip production 
cross section. The gg components of the Xcj cross sections, all singlets, are 



^{99 
^{99 
^{99 



XcO) 



Xcl, 



Xc2) 



-B26,- 



lXci/3 , 



--D3t7a: 
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15 



+ 



h{z) 



4^2 In zfi{z) 

1 + Z)5(1-Z)4 ' 3(1 + ^)4(1-^)3 



5 B{99 Xco) 



(19) 
,(20) 
(21) 
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where 



(z) = + 9z'^ + 26z^ + 28z^ + 17 z^ + 7z^ - AOz^ - 4^ - 4 , 
f^^z) = + 39z^ + Ubz"^ + 251z^ + 119z^ - 153z^ - 17z^ - U7z^ - 8z + 10 



The Xci also has a singlet contribution from gq scattering, 



^{gci Xci) = 777 ^s^x 



8 o . {ore Pi)) 



1 " " 



—z'^ In z + 



4z-^ -9z + 5 



(22) 



The qq contributions to XcJ production are all color octets, 



1 fi 

Hqq^Xco) = -B,5,{Or('S,)) , (23) 
^ Xcj) = (2 J + l)a(gg Xco) • (24) 

The sjTumetry relations have been used to obtain the gg contribution to Xc2 produc- 
tion, Eq. (pT[), and the qq contributions to Xa and Xc2 in Eq. (p^. The relevant 
nonperturbative matrix elements for Xcj production are 

(g*! i^o)) ^ ^^^3 ^ /(>^xco N\ _ Q QQ32 Qgv^ _ (25) 
In Ref. the singlet matrix elements in Eqs. (|l^) and (|25|) were calculated from 



the quarkonium wavefunctions at the origin. The octet matrix elements were fit to 
Tevatron production data and Ag was obtained from a fit to total cross sections data 
at fixed-target energies. Note that in NRQCD, three parameters are needed to fix 
the ip' production cross section while eight are needed for the total cross section. 
Only one parameter for each state is needed in the GEM, a considerable reduction. 

The total ip forward xp distributions^ at 800 GeV and 120 GeV, Eq. (|T2[), are 
shown in Fig. |^(a) and (c) respectively. Since the Xco branching ratio to is less 
than 1%, its contribution is virtually negligible. However, ^ 27% of the produced 
Xci states decay to ip, including the gq scattering contribution, Eq. (0), shown in 
the dotted curves. At 800 GeV this component is only a factor of 2-3 less than the 
total qq contribution to the full ip cross section. The gg contribution from Xci decays, 
Eq. (^0|) , and the smaller Xc2 decay contribution, 14%, provide most of the singlet 
component of total ip production. Interestingly, when the Xcj decays are included, 
the octet contribution to the total ip production cross section is 60%, close to the ip'. 
At 120 GeV, the percentage of octet production is larger for both the total ip and 
ip', 67% and 78% respectively. The ip' distributions at 800 GeV and 120 GeV are 
given in Fig. ^(b) and (d), including the individual contributions from gg fusion and 
qq annihilation. Here the qq contribution is largest at xp ~ 0.7 at 800 GeV and ~ 0.6 
at 120 GeV, similar to the GEM. However, for the total ip distribution, the qq 
contribution is significantly smaller than in the GEM, resulting in a narrower ip xp 
distribution in NRQGD. The NRQGD cross section is a factor of 3 smaller than the 



"'Note that, as in the CEM, the pp xp distributions are symmetric around xp — 0. 
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Figure 3: The V xp distributions at (a) 800 GeV and (c) 120 GeV in NRQCD. The 

contributions from gg fusion (dashed), qq annihilation (dot-dashed), gq scattering 
(dotted) and the total (solid) are given. The corresponding ip' distributions arc given 
in (b) and (d). The curves show gg fusion (dashed), qq annihilation (dot-dashed), 
and the total (solid). 
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CEM cross section at xp ~ 0.9 and a factor of 1.5 smaller at xp ~ 0.5 where the two 
are equal at Xi? = 0. The qq component does not dominate the total ip distribution 
until xp ~ 0.8 at 800 GeV and ~ 0.7 at 120 GeV. These relative differences can 
influence the strength of nuclear effects that depend on the nuclear quark and gluon 
distributions. 



3 Drell-Yan Production 

Lepton pairs are produced by the Drell-Yan process, qq annihilation into a virtual 



photon at leading order, qq ^ ^* ^ l^l |^. The partonic cross section for Drell- 
Yan production is 

-el5{s-M^) (26) 



dM 9M 

where s = X1X2S. To obtain the hadroproduction cross section as a function of 
pair mass, M, and xp, we must fold the partonic cross section with the quark and 
antiquark densities evaluated at M, here taken to be in the range 4 < M < 9 GeV, 
between the ip and T family regions. Then 



/ dxidx2 5{xiX2S — M^) 5{xp — xi + X2] 
Jo 



dxpdM ~ 9M Jo 

X Ee;[/n^i'^')//(^2,M2) + /|(xi,M2)/^^(x2,M2)] . (27) 



After integrating the delta functions, the LO cross section, including the isospin of 
the target nucleus, is 

+ /|(i„i, Af2)(z^/P(i„2, ) + ru/,»(i02, M^))] , (28) 

where za = Z/A and ua = N/A are, respectively, the fractions of protons and 
neutrons in the target nucleus. 

When this leading order cross section is compared to data, it falls short by an ap- 
proximately constant factor, known as the K factor. Experimentally, it is ~ 1.7 — 2.5, 
depending on the energy, mass range, and parton distribution functions. At NLO, the 
Compton and annihilation processes, qg —>■ q'-f* and qq —>■ (77* respectively, contribute 
in addition to virtual corrections to the LO cross section, resulting in a theoretical K 
factor — the ratio of the NLO to the LO cross sections — of approximately 1.4 — 2, 
somewhat less than that obtained by comparison to the data [^]. This theoretical 
K factor serves the same purpose as the adjustment of between the LO and NLO 
calculations in the CEM, discussed earlier. 

In Fig. ^(a) and (c) the Drell-Yan xp distribution is shown for masses between 4 
and 9 GeV at 800 GeV and 120 GeV. The 120 GeV xp distribution does not extend 
over all xp since the phase space for high xp and high mass pairs becomes limited. 



11 




0.00 0.25 0.50 0.75 4 6 8 



Xf M (GeV) 

Figure 4: The Drell-Yan xp distributions for 4 < M < 9 GeV at (a) 800 GcV and 
(c) 120 GeV. The Drell-Yan mass distributions, integrated over xp are shown in (b) 
and (d) at 800 and 120 GeV respectively. The leading order results are given by the 
solid curves, the next-to-leading order results are shown in the dashed curves. 
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Because the quark distributions have a harder x dependence than the gluon, the 
Drell-Yan xp distribution is broader than the ip xp distributions shown in Figs. ^ and 
^. For comparison, both the LO and the NLO distributions are shown. There is some 
dependence of the K factor onxp- At 800 GeV and ~ 0, the theoretical K factor is 
1.4, increasing to 2.1 at xp = 0.9. The calculated K factor is 20% larger at 120 GeV. 
The change with xp reflects the increasing importance of the Compton process with 
increasing xp, corresponding to an increase in the gluon density at low X2- The Drell- 
Yan mass distribution, integrated over the corresponding xp ranges given in Fig. ^(a) 
and (c), is shown in Fig. ^(b) and (d) at 800 and 120 GeV respectively. Since xp is 
integrated over, the lower xp values are the most important for the determination of 
the K factor as a function of mass. The theoretical K factor is ~ 1.4 at 800 GeV and 
1.8 at 120 GeV and does not change more than a few percent with mass. 

Since the Drell-Yan mechanism produces lepton pairs which only interact elec- 
troweakly, the A dependence is expected to be weak because no final-state interac- 
tions affect the lepton pair. However, initial-state interactions such as shadowing and 
energy loss may influence the A dependence, as we discuss in Sections 6 and 7. 

4 Nuclear Absorption in pA Interactions 

The cc pair may interact with nucleons and be dissociated or absorbed before it can 
escape the target. The effect of nuclear absorption alone on the ip production cross 
section in pA collisions may be expressed as 

cTpA = cjpN I (fhTfih) , (29) 

where h is the impact parameter and T'^{h) is the effective nuclear profile function, 

/oo 
dzpA{h,z)S--''' . (30) 

The probability for the cc pair to avoid nuclear absorption and form a if), called the 
nuclear absorption survival probability, S^^^, is 

^abs = exp |- dz'pAib, z')cr,U^' - ^)| (31) 

where (Tabs is the charmonium (or ccg |^) nucleon absorption cross section. The 
nuclear density profile is T^ib) = J^^dzpA^b, z) so that T'^ib) = TAib) when S^'^^ = 
1. Nuclear charge density distributions from data are used for pA [^. Note that 
expanding S'^^^, integrating Eq. (|29D , and reexponentiating the results assuming A is 
large leads to Eq. (|I|) with a = 1 — 9(Tabs/(16vrrQ). 

We consider three different models of nuclear absorption: either all quarkonium 
states are produced as color octets or color singlets or as a combination of octet 
and singlet states. When pure octet or pure singlet absorption is considered, ip and 
ip' production are calculated in the GEM. When a combination of octet and singlet 
absorption is assumed, NRQGD is used to obtain the correct balance between octet 
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and singlet production. Both the CEM and NRQCD model parameters are tuned to 
fit pp production. In this section, we only give examples of a range of cross sections 
for each absorption model. The actual values of (Tabs are set in Section 9 after initial 
state effects have also been included. 

In our considerations of absorption, we include the ~ 30% contribution from 
Xcj decays |27| and the ~ 12% contribution from ip' decays ||2l[] decays. Then the 



total ip survival probability, including indirect production, is 

= 0.5SS^% + 0.35 + 0.125;^^^ . (32) 

The ip' itself is only produced directly since other, more massive, charmonium reso- 
nances decay to DD pairs. 

The first case, pure octet production, assumes that all charmonium states are 
initially produced as \ccg) states with the same absorption cross sections, leading 
to (Jabs = cr^N — ^i)'N — ^xcjN- Therefore the survival probabilities, Eq. (pTD, are 
identical for all states. Thus the feeddown contributions to the ip iii Eq. ( P^ do 
not affect the absorption. After ~ 0.3 fm/c, the remaining \ccg) states are expected 
to hadronize. Since the absorption cross section is established at the production of 
the state, the octet cross section is independent of the position z and thus xp and 
projectile energy. In this model we treat absorption as if only the \ccg) interacts with 
nucleons, not the final charmonium states. In Fig. §(a), a is given for several values 
of the \ccg) cross section: dabs = 1, 3, 5, and 7 mb corresponding to a; = 0.98, 0.95, 
0.92, and 0.90 respectively. It is obvious that octet production alone will not modify 
the shape of as a function of xp- 

We now discuss absorption when all cc pairs are produced as color singlets. If 
the cc pair is produced as a color singlet, it is initially small with a spatial extent on 
the order of its production time, r oc m~^, ignored in the calculation. The proper 
time required for the formation of the final charmonium bound state obtained from 
potential models [^, ~ 1 — 2 fm, is considerably longer. The cc-N absorption 
cross section may be expected to grow as a function of proper time until r^- when 
it saturates at the asymptotic value cr^.^- We simulate the growth of the absorption 
cross section by ^ 



(Tshsiz - z) = { \T^J . (33) 

7^.^ otherwise 

The exponent k determines the increase of (Tabs during hadronization of the cc pair. 
If (Jabs is proportional to the geometric cross section, then we expect k ~ 2. (See 
also for predictions of a^^ if k = 1.) The proper time r is related to the path 
length traversed by the cc pair through nuclear matter, r = [z' — z^j'yv. The 7 factor 
introduces xp and energy dependencies in the growth of the cross section. Depending 
on the initial energy of the projectile and the size of the target, the cc pair may form 
a '(/' inside or outside of the target. 

In Fig. |^(b), examples of a are given for direct ■?/' and ip' absorption at 800 GeV. 
The solid and dashed curves assume (j^^ = 5 and 10 mb respectively while the dot- 
dashed and dotted curves are for cr^/^v = 1^ 20 mb. In our calculations, we assume 
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Figure 5: The A dependence of nuclear absorption models is given in (a), (b) and 
(c) and the comover A dependence is shown in (d). In (a), octet cross sections of 1 
mb (solid), 3 mb (dashed), 5 mb (dot-dashed) and 7 mb (dotted) are shown. Singlet 
absorption is shown in (b) for ip with cr^jy = 5 mb (solid) and 10 mb (dashed) as 
well as ip' with a^,pf — 15 mb (dot-dashed) and 20 mb (dotted). A combination of 
octet and singlet production is assumed in (c). The curves represent: ip absorption 



with fx-t/t 



1 mb and a. 



singlet 
abs 



1 mb (solid) and cr^i^J^^ — 3 mb and a 



(dot-dashed); ^p' absorption with a°^l^ 



cr°^s^* = 3 mb and a. 



singlet 
abs 



1 mb and (T^J^g^'*^* 



singlet 
abs 



5 mb 



3.7 mb (dashed) and 



19 mb (dotted). In (d), comover interactions are shown 



for a^co = 0.67 mb (solid) and cr^'co = 3.7(T^co (dot-dashed). 
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that the asymptotic absorption cross sections scale in proportion to the squares of the 
charmonium radii cr^ij^ ~ 3.7(7^^ and o"^^_,Ar ~ 2.4(Tf,,Ar. Thus each contribution 



to Eq. ( |32D has a different A dependence. The ip and ip' formation times are different, 
= 0.92 fm and t^> = 1.5 fm [0]. The ip and ip' results at xp < in Fig. |^(b) reflect 
the differences in formation times as well as the gamma shift due to their masses. At 
800 GeV, hj xp = the final-state meson is produced outside the target so that 
a ^ 1 for xp > 0. Therefore the A dependence of color singlet production is virtually 
independent of ct^.at for x^? > at 800 GeV. At 120 GeV, both states can be produced 
inside the target at xp > and influence the A dependence at forward xp as well. 

More realistically, ip production is a combination of octet and singlet states, as in 
NRQCD. The ratio of octet to singlet production is energy and xp dependent |^ so 
that the relative absorption of each state depends on xp since the octet and singlet 
absorption cross sections are expected to be different [|^. Because the ip' is directly 
produced, the xp dependence of absorption is straightforward. 



da, 



pA 



dxi 



J i/)',oct 

dxp 



d%T 



oflf (oct) 



J V'. sing 

dxp 



d%T 



cfT (sing) 



ib), 



(34) 



where 0"°^* and cr^™g^ replace a^hs in Eq- (POD- The ip distribution is more compli- 
cated since we must account for the fact that the octet absorption cross section is 
independent of the charmonium state while the singlet cross sections are not. In the 

-lefT (oct) , 
A \ 

absorption is different for each individual state. Then, 



octet case, the same can be applied to all states feeding the ■?/' while singlet 



7 ih, tot 

d<^pA 



dx. 



J i/), dir,oct 2 ^^-XcJ. oct 

CLX p J dx p 



J=0 
J ^'4'', oct 

dxp 

dir,sing 



^2^ycfr(oct) 



(b) 



+ d'b 



pp 



dxi 



yj.dir,eff(sing)^^^^^^^^^^ 



J=0 



+ B{iIj' ^ iIjX)- 



pp 



1-0', eff (sing) 



dxf 



(b) 



(35) 



pp 



dxi 



_rpXc.J,C^ '■'^™^-'(6) 



In Ref. EH], the singlet cross section was assumed to be negligible. Therefore, a. 



oct 
abs 



11 mb was needed to produce an effective a equivalent to the assumption of pure 
octet production with (Tabs = 7.3 mb obtained in |TT[]. We use Eq. ( PBj ) and cr^^ ^ 0. 



Also, in Ref. the authors only calculated the x^-integrated cross sections. Here 
we use Eqs. (|^) and (|3l) with the full xp dependence to calculate a{xp) for ip and 
'0' absorption, shown in Fig. |^(c) at 800 GeV. The differences between the ip and ip' 
results arise from their distinct xp dependencies in the NRQCD model, predominantly 
from the XcJ contributions to ip production. 

We have illustrated several different combinations of ip' and ip absorption cross 
sections in Fig. §(c). We choose cr^?^* = cr'^}f^ as in pure octet production. Fig. §(a), 
and ct^^jv^'^* ~ 3.7cr^^Jy as in the pure singlet case shown in Fig. ^(b). The differences 
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in for and ip' are small, especially between the solid and dashed curves with 

"^abs*^* = 1 mb. The only obvious differences are aX xp < when the ip' singlet 
contribution is larger and the ip' is still produced inside the target. However, at large 
\xf\i qq annihilation, an octet contribution to begins to become more important, 
causing the change in slope of a{xF) here. This effect is not as strong for the i/j because 
the qq contribution does not overtake the gg until larger xp. It is interesting to note 
that the effective a in the solid curve (Cabs*^* = 1 and o"^™^^^'^* = 1 mb) is similar to 
the 

'''abs ~ 1 result shown in Fig. ^(a) for pure octet absorption. Assuming a 3 mb 
octet absorption cross section results in a similar effective a at a; i? < in Figs. |^(a) 
and (c) — compare the dashed curve in (a) with the dot-dashed curve in (c). However, 
at forward xp, the effective a in (c) is larger since the color singlet components 
escape without absorption, at least until the growing qq contribution causes the octet 
mechanism to dominate ip absorption once again at xp > 0.7. The effect of singlet 
absorption would be even weaker at 120 GeV because the octet contributions make 
up a larger fraction of the production cross section at this energy. 

The results with the absorption cross sections shown in Fig. |^ are only examples 
of the magnitude of the effects. It is clear, both from the data in Fig. |^ and from 
the initial state effects discussed in the following sections, that the model absorption 
cross sections must be smaller than those used previously when no initial state effects 



were included 10, 11 



5 Hadronic Comovers in pA Interactions 

Comoving secondaries, formed after tq ~ 1 — 2 fm, may also scatter with the cc pair 
or the ip. Because ~ tq, the final-state charmonium is assumed to interact with the 
comovers. A spectator hadron moving with a velocity close to that of the charmonium 
state enhances the dissociation probability. 

The A dependence of ip production due to comovers alone is determined from 

<JhA = <JhN J d%S^\h) , (36) 

where the total probability that the ip survives its interactions with comovers is 

= 0.58^-,i, + 0.35-^ + 0.125- . (37) 



The direct ^-comover survival probability is [|1^ 

-^^dirl^) ~ exp |- j dr {a^covMr, 6)| . (38) 

The other survival probabilities for comover interactions with charmonium states 
are similar. The parameters are the charmonium-comover absorption cross sections, 
the velocity of the relative to the comovers, v ~ 0.6, and n{T,b), the density of 
comovers at time r and impact parameter b. We take a^co = 0.67 mb from a study of 
ip suppression in nucleus-nucleus data [|T0| with a^'co ~ 3.7(j^co and Cxcjco ~ 2.4cr^co 



44], assuming that the asymptotic charmonium states interact with the comovers. 
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Integrating Eq. ([ 

hadron rapidity density. 



over T and relating the initial density of the system to the final 



noTo 



-^{dN/dy) H, one finds p 



dTn{T, b) 



In 



ri\dN 

— -i-(^hNTA{b) 
To J dy 



(39) 



where the effective proper lifetime tj over which the comovers interact with the ip is 
'Ti ~ ^p/v and Tp ~ 0.8 fm. This scaling assumption is rather strong for pA since no 
collective motion is expected. However, perhaps within the interaction tube carved 
out by the incident proton, scaling may hold. There is evidence of rapidity scaling 



for produced particles in the central region of pp collisions |57 . 

The rapidity density grows with center of mass energy ||4^ . The shape of the pro- 
duced particle rapidity density with inclusive ip production is unknown. We assume 
that the multiplicity slope is the same on both sides of midrapidity 



dN 
dy 



dN, 
dy 



\y=o - ay , 



(40) 



47| 



where we take dN/dy\y=Q = 1.07 in pp interactions at 800 GeV and a = 0.108 
The CO mover density is depleted at forward rapidities but enhanced close to the target. 

Since the transverse area over which the comovers are produced, ttR^, is approxi- 
mately equal to ChN, the direct ip survival probability in hA collisions may be recast 
as 



oco 

tp, dir 



(41) 



The similarity between Eqs. (^91) and (^iD suggests that ■j/'-comover interactions do not 
introduce any unusual A dependence. Thus comover contributions to pA interactions, 
while small, are difficult to rule out entirely. 



6 Nuclear Shadowing 

Measurements of the nuclear charged parton distributions by deep-inelastic scattering 
off both a large nuclear target and a deuterium target show that the ratio = 
F2 / F2 has a characteristic shape as a function of x |jl2|. The region below x ~ 0.1 
is referred to as the shadowing region and the range 0.3 < x < 0.7 is known as the 
EMC region. In both regions, the parton density is depleted in the heavy nucleus 
relative to deuterium, i.e. < 1. At very low x, a; ~ 0.001, appears to saturate 



i8[] . Between the shadowing and EMC regions, an enhancement, antishadowing, is 
seen where Rf2 > 1- There is also an enhancement as x — > 1, assumed to be due 
to nucleonic Fermi motion. The general behavior of Rf2 as a function of x is often 
referred to as shadowing. Although this behavior is not well understood for all x, the 
shadowing effect can be modeled by an A dependent fit to the nuclear deep-inelastic 
scattering data. 

We have assumed that the nuclear parton distributions factorize into the nucleon 
parton distributions, independent of A, and a shadowing function that parameterizes 
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the modifications of the nucleoli parton densities in the nucleus, dependent on A, x, 
and Q^: 

ff{x,Q\A) = S\A,x,Q')fr{x,Q') . 



While the location of the parton in the target could influence 5* the impact pa- 
rameter is difficult to resolve in pA collisions. We use three different parameterizations 
of the shadowing function, S\A,x,Q^). 

The first parameterization is a fit to nuclear deep-inelastic scattering data which 
does not differentiate between quark, antiquark, and gluon modifications and does 
not include evolution in Q^. Therefore it is not designed to conserve baryon number 
or momentum. We define = Si{A,x) |Q with 



Si{A,x) 



l + 0.0134(l/x-l/xsh) 



1 _ ■ \ 0-321 



(42) 



R f ( ^ _ ^ ) Xfermi < X <1 



where R^ = aemc-^emca^sh, Rf = aemc-^cmca^fermi, Kmc = 0.525(1-^-^/^-1.145^-2/3 + 

0.93A-^ + 0.88^-^/3 - 0.59A-5/3), and Oemc = 1 + ^emca^emc- The fit fixes the x values 
at the boundaries of the x regions, Xgh = 0.15, Xemc = 0.275, and Xfermi = 0.742. Thus, 
the nuclear parton distributions are modified so that 

ffix,Q') = 5i(Ax)/f(x,Q2). (43) 

The parameterization is available for all A and is designed so that 5*1 = 1 when 
A = 1. Figure ^(a) shows the parameterization for A = 184 and A = 9. Note that the 
antishadowing region is rather narrow and saturation appears at x < 10-^. Figures 
|^(b) and (c) give a^xp) for ip and Drell-Yan production respectively. Since this 
parameterization affects all partons equally, the results are independent of the chosen 
parton distribution function. They are also virtually independent of the charmonium 
production mechanism although there is a slight model dependence because the CEM 
involves an integral over Imd \fs < x < 2mo/y/s, Eq. (P), while in NRQCD xi and 
X2 are either both fixed, as in Eqs. (pTSD, ([T^), (plf), (p^), and (0) or xi is fixed by 
the delta function in Eq. (10) while Am'^/xiS < X2 < 1. 

The second parameterization, 5*2 (A, x, Q^), modifies the valence quark, sea quark 
and gluon distributions separately and also includes evolution beginning at 



Q = Qo = 2 GeV and continuing up to Q = 10 GeV. It is based on a fit to the data 
using the Duke-Owens parton densities [^. In this case, the nuclear parton densities 
are modified so that 

f^{x,Q') = S^{A,x,Q')f^{x,Q') , (44) 
fi{x,Q') = SUA,x,Q')fs{x,Q') (45) 
f^{x,Q') = S^{A,x,Q')fai^,Q') , (46) 

where fv = uv + dy is the valence quark density and fs = 2(u + d + s)is the total 
sea quark density. It is assumed that S2 and are the same for all valence and sea 
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Figure 6: (a) The shadowing parameterization for W (sohd) and Be (dashed) 
targets as a function of x. The resulting A dependence for (b) ip production in 
the CEM and (c) Drell-Yan production is given at 800 GeV (sohd) and 120 GeV 
(dashed). The NRQCD ijj results are shown in (b) at 800 GeV (dot-dashed) and 120 
GeV (dotted). 
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quarks, consistent with the symmetric sea of the Duke-Owens parton distributions. 
These modifications conserve baryon number, Jq dx fy{x,Q'^) = Jq dx fy{x,Q'^), and 
the parton momentum sum, X]p ^fpi^iQ"^) = J2p lo dx xfp{x,Q'^) where P = 
V, S, and G, at all Q^. Using parton densities other than Duke-Owens may lead to 
small deviations in the conservation rules. 

The parameterization is only available for A = 32 and 200. It is thus applied only 
to the tungsten target and the beryllium densities are left unmodified. Figure |^(a) 
shows the ratios 5^, and 5*^ at Q = Qo and Q = 10 GeV. At Qo the sea quarks 
are shadowed more strongly at low x than the gluons. Both the valence quarks and 
gluons are antishadowed while the sea quarks are not. The effects of evolution are 
weakest for the valence quarks and strongest for the gluons. Figures ^(b) and (c) 
show a{xF) for ip and Drell-Yan production respectively. The ip results in the CEM 
are given for the MRST LO pSf distributions and in the NRQCD approach with the 
CTEQ 3L parton distributions. The main differences in the production models 
appear at negative xp, corresponding to the EMC dip at large x and appears because 
of the evolution of the gluon distributions at large x. The two calculations evolve 
differently because rric = 1.2 GeV in the CEM and 1.5 GeV in NRQCD. The larger 
scale causes a smaller EMC dip in the shadowing ratio for the NRQCD calculation. 
The differences between production models at large xp are due in part to the gq 
scattering contribution. Since this component is virtually negligible at 120 GeV, the 
model dependence is then small for 5*2. Choosing other parton distribution functions 
for CEM ip and Drell-Yan production results in very similar ratios as for MRST LO. 

A more recent shadowing parameterization, Sl{A,x,Q'^), based on the GRV LO 
parton distributions ||5^, is now available []54| , The initial scale was chosen to 



equal the charm quark mass in the GRV LO distributions, Q = Qo = 1.5 GeV. At this 
scale all sea quark ratios are assumed to be equal, as are both the valence ratios. The 
parameters are constrained by nuclear deep-inelastic scattering and Drell-Yan data. 
The gluon ratio is then fixed by the momentum sum rule as well as ip electroproduction 
data. Above Qq, the individual quark and gluon distributions are evolved separately. 
The gluon distribution has a larger antishadowing peak in this parameterization while 
the sea quarks are shadowed in the same region, a significant difference from 5*2. The 
Drell-Yan data on the violation of the Gottfried sum rule is taken to account, thus 
5*3 7^ 5*1 above Qq. Evolution is taken up to Q = 100 GeV and the parameterization 
is generalized to all A, both improvements over 5*2. Again however, using other parton 
densities besides GRV LO could lead to small deviations from the conservation rules. 

In Fig. 1(a) and (b), we show ratios for the uy, u and g densities at Q = Qo and 10 
GeV for W and Be respectively. Figures ||(c) and (d) give the corresponding a{xp) for 
ijj and Drell-Yan production. The tp results in the CEM are given for the MRST LO 
distributions. Since gluon shadowing is not as strong as in the 5*2 parameterization 
at low X, the effective a is larger at large xp than in Fig. 0(b). We have checked the 
CEM results with other parton distributions and found that the differences between 
the parton distributions are also more pronounced at large xp since the individual 
quark and antiquark distributions evolve separately while with 5*2, the valence and 
sea quarks respectively were considered together. In NRQCD, the effective a with 
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Figure 7: (a) The 5*2 shadowing parameterization for A = 200 as a function of x. The 
valence ratios, Ry, are given by the sohd curves, the sea quark ratios, Rs, by the 
dashed curves and the gluon ratios, Rq, are given by the dot-dashed curves. At small 
X, the lower curves arc at Q = 2 GeV and the upper are at Q = 10 GcV. In (b) the 
ip A dependence is illustrated for the GEM with MRST LO distributions at 800 GeV 
(sohd) and 120 GeV (dot-dashed). The NRQGD results with the GTEQ 3L densities 
are also shown at 800 GeV (dashed) and 120 GeV (dotted). In (c) the Drell-Yan A 
dependence is given for the MRST LO distributions at 800 GeV (sohd) and 120 GeV 
(dashed). 
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Figure 8: The ^'3 shadowing parameterization as a function of x for (a) W and (b) Be 

targets. The valence up ratios, Ruy, are given by the sohd curves, the u ratios, Ru, 
by the dashed curves and the gluon ratios, Rq, are given by the dot-dashed curves. 
At small X, the lower curves are at Q = 1-5 GeV while the upper curves are at Q = 10 
GeV. The resulting A dependence for (c) ip and (d) Drell-Yan production is given. 
In (c) the t/j A dependence is illustrated for the CEM with MRST LO distributions 
at 800 GeV (solid) and 120 GeV (dot-dashed). The NRQCD results with the CTEQ 
3L densities are also shown at 800 GeV (dashed) and 120 GeV (dotted). In (d) the 
Drell-Yan A dependence is given for the MRST LO distributions at 800 GeV (solid) 
and 120 GeV (dashed). 
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the CTEQ 3L densities at 800 GeV is similar to the GEM results except at larger xp 
due to the gq scattering contribution, as with the 5*2 parameterization. However, at 
120 GeV, the model dependence is more pronounced than with 5*2, due to the larger 
relative importance of qq annihilation in the GEM than in NRQGD. This difference 
is again less important with the 5*2 parameterization because it does not distinguish 
between the individual quark and antiquark distributions. At negative xp there is no 
difference due to evolution at the EMG dip because the ratios at large x in Figs. |^ 
(a) and (b) are essentially independent of Q^. The Drell-Yan results are only shown 
for the MRST LO distributions. The reduced antiquark shadowing at low x results 
in a larger a than with the 5*2 parameterization. 

To summarize, we note that the shape of a{xp) is fixed by each parameterization. 
It is clear from the results in Figs, j^-]^ that shadowing alone is insufficient to describe 
the preliminary E866 ip data as a function oi Xp- This fact has been known for some 
time since the NA3 ||^ and E772 |p A dependence was similar as a function of 
xp but not as a function of X2 as would be expected if the nuclear dependence was 
dominated by shadowing. 



7 Effects of Energy Loss 

Partons are expected to lose energy when traversing matter. This effect has been 



discussed primarily in the context of jet quenching [58[] . Since the projectile 
parton is typically expected to feel the effects of energy loss, the scaling of the A 
dependence at different energies with xp oi xi suggested that energy loss could be 
the cause. We will introduce three models of energy loss that have been applied 
earlier to ip production and discuss their influence in the context of the E866 data. 

7.1 Initial State Loss 

Initial state energy loss, as studied by Gavin and Milana [|13| and subsequently devel- 
oped by Brodsky and Hoyer |T^, takes a multiple scattering approach that essentially 



depletes the projectile parton momentum fraction, xi, as the parton moves through 
the nucleus. Both the quarks and gluons can scatter elastically and lose energy before 
the hard scattering. This loss produces a similar effect for Drell-Yan and ip produc- 
tion. The motivation for this model stemmed from the fact that the A dependence 
of ip production at 200 and 800 GeV seemed to scale with xp (or xi) and not X2 
P, 1^. The projectile parton momentum fraction involved in the hard scattering is 
then x'l = Xi — Axi where xi is the original projectile parton momentum fraction 
when the parton first entered the target and Axi represents the loss in xi due to 
multiple scatterings. Thus the shifted value, x[, enters the partonic cross sections 
but the parton distributions must be evaluated at the initial Xi. An additional delta 
function is added to Eqs. (^) and (|TD|) with the corresponding integral over x'l so that 
Eq. (H) becomes 



da""" 
dxpdrn? 



1 

- / dx[dxidx2S{x[ — Xi + Axi) 
s Jo 
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X S{xf — x[ + X2) S{x[x2S — m^) Hab{xi,x[, X2; m?) (47) 

while Eq. (|l^) is then 



dxp 



j dx[dxidx2 6{x'i — xi + Axi) 
X 5{XF ~ x[ + X2)ff{xu fi')ff{x2, fJ^'Mtj ^ C) . (45 



We first discuss the model by Gavin and Milana |13] and then the modifications 



suggested by Brodsky and Hoyer [T^ with later refinements by Baier et al. |^ . 



The first model of initial-state energy loss apphed to J /ip production was proposed 
by Gavin and Milana referred to as GM hereafter. In their model, they assumed 
that 

Axi = e,XiAi/='(^^j'" (49) 

with n = 1. We do not include the dependence in our calculations so that here 
n = 0. The energy loss depends on the parton identity in this formulation. The initial 

Xi is 

where i = q or g with = 0.00412 and tg = 9eg/4 due to the difference in the 
color factors. When n = 0, Eq. (p9D corresponds to —dE/dz\q ~ 1.5 GeV/fm and 
—dE/dz\g ~ 3.4 GeV/fm |T^. In our calculations, we assume only initial state 
elastic scattering of the quarks and gluons. Final state effects on the if) included in 



Ref. |13[ are left out here under the assumption that final-state absorption provides 
a compensatory effect. 

In Fig. P we show the results for this mechanism alone on the A dependence of 
and Drell-Yan production. The A dependence is weak at negative xp where xi is 
already small so that further reduction does not significantly change the quark and 
gluon distributions. This is true even for parton distributions that increase as xj"" 
when Xi is small and a = 0.3 — 0.5. As xp increases, Xi grows larger and if the parton 
densities behave as ~ (1 — Xi)"^ as Xi — 1, a slight decrease in xi is magnified. The 
effect should be stronger for than Drell-Yan production because ~ 5 > Ug^- ~ 3 
in simple spectator counting models ||59| and the valence quark distributions are 
most important for Drell-Yan production at large xp (and Xi). The choice of parton 
densities does not change the shape of a^xp). The energy dependence is also rather 
weak. A comparison of Figs. |^(a) and (b) shows that the behavior of a{xp) does not 
depend strongly on the production model although there is evidence that the effect 
begins to be nonnegligible at a lower Xi;' in the NRQCD approach. 



Later, Brodsky and Hoyer [Q, BH, argued that the energy loss in the Gavin and 



Milana model was too large because there is not enough time after the initial QCD 
bremsstrahlung for the color field of the parton to be regenerated. Therefore, the 
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Figure 9: The A dependence assuming GM loss for ip production in (a) the CEM 
and in (b) NRQCD and (c) Drell-Yan production. In (a) the ijj A dependence is 
illustrated for the MRST LO distributions at 800 GeV (solid) and 120 GeV (dashed). 
The NRQCD results with CTEQ 3L are shown in (b) at 800 GeV (sohd) and 120 GeV 
(dashed). In (c) the Drell-Yan A dependence is given for the MRST LO distributions 
at 800 GeV (sohd) and 120 GeV (dashed). 
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subsequent interactions of the parton in the target do not lead to a large increase in 
energy loss ||6^. From the uncertainty principle they deduced that the loss should be 



independent of parton type and the change in Axi should be bound so that 

A.. < <^ (51) 

where La is the path length through the medium and {k\) is the average transverse 
momentum of gluons radiated by the incoming parton. \i E = xis/2mp and La ~ 
Ra oc then 

Axi < —A'l^ (52) 

XiS 

where n oc mp{k'j_). The average radiative loss is thus expected to be —dE/dz ~ 0.25 
GeV/fm with another 0.25 GeV/fm loss expected to arise from elastic scattering. In 



this case, when Axi oc c/xi, Xi = 0.5{xi + \J{x'iY + 4c). The xp dependence of a 

when c = kA^/^/s, referred to henceforth as "original BH loss", is given by the dotted 
and dot-dash-dashed curves in Fig. ^ for ip and Drell-Yan production. 

Subsequently, the bound on dE/dz was refined through the work of Baier et al. 
Il58| , |6l|| where they determined 

with {p\w) the characteristic squared transverse momentum of the parton|. The 
value of the radiative loss is independent of the details of the scattering process as 
long as La is large. In this description, Axi is then 

Axi = ^-^^La{pI^) (54) 
2 XiS 

where the average transverse momentum is proportional to A^/'^ [Q. Since 

{p\w) ^ A^^^^ Axi oc A^f'^ in Eq. (0) rather than A^/'^ as postulated by Brodsky 
and Hoyer |]T3[, Eq. (|52D, because they assumed that {k\) was independent of A. 



Two estimates of {p\\y) were provided in Ref. [Q. The larger value, used as an 
upper limit, comes from a single nuclear rescattering of photoproduced dijets ||63[| , 



(Pin-) = 'r^a-Ai^.^'^^^i^S^i^i^ . (55) 

They obtained 0.05 < Alqs < 0.1 GeV^ by assuming that dijet production is domi- 
nated either by quarks or gluons using the measured pt broadening as a function of 
A. With the lower bound on Alqs, 

(p^^) ~ 0.658 as A^'^ GeV^ . (56) 



^ After the inclusion of other diagrams suggested by Zakharo v [|62[ , Baier et al. concluded that 
the loss derived in Ref. |Q was, in fact, a factor of two larger |6l[|. This difference is reflected in 
Eq. dH). 
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Figure 10: The A dependence assuming BH loss for ip in (a) the CEM and in (b) 
NRQCD and (c) Drell-Yan production. The MRST LO distributions are used for 
CEM and Drell-Yan production while CTEQ 3L densities are used with NRQCD 
ip production. The maximum BH loss, Eq. (p6D , is shown in the solid curves at 800 
GeV and in the dot-dashed curves at 120 GeV. The minimum BH loss, Eqs. (^) and 
(pSD, is shown in the dashed curves at 800 GeV and in the dot-dot-dot-dashed curves 
at 120 GeV. The original BH loss, Eq. (|^), is shown in the dotted curves at 800 GeV 
and in the dot-dash-dash-dashed curves at 120 GeV. 
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Since the initial states could not be explicitly identified, we assume that {p]_w) is 
identical for quarks and gluons. Then when ~ 0.3 and A = 184, we find —dE/dz ~ 
1.28 GeV/fm with Eq. (|56D. We refer to this as "maximum BH loss" in the remainder 
of the discussion even though dE/dz is actually smaller than the proposed GM loss. 
(The difference in the subsequent shapes oia{xF) lies in the form of Azi.) The second 
estimate depends on the nucleon gluon distribution and contains explicit color factors 
so that 

{pIw), = ^^PAxG{x,Q')LA^0.07a,A'/-'GeV' (57) 

{Plw)a = ^(plvK).^0.15a,AV3GeV2 (58) 

where xG{x) ~ 1 — 2 for the xi range of E866. This lower estimate is referred to 
subsequently as "minimum BH loss". Now when as ~ 0.3 and A = 184, —dEq/dz ~ 
0.12 GeV/fm and -dEg/dz ~ 0.28 GeV/fm. 



The resulting xp dependence of a is shown in Fig. |T0| at 800 and 120 GeV. At 
negative xp-, Xi can be considerably larger than x[, up to an order of magnitude as 
Xi? — >■ — 1 at 800 GeV. The difference in shapes at negative xp between the GEM and 
NRQGD arise from the relative importance of qq annihilation and gg fusion (as well 
as qg scattering in NRQGD). Even though the energy loss is the same for quarks and 
gluons with the original and maximum BH loss estimates, the relative change is larger 
for the gluon than the sea quark distributions when xi is small. At large negative 
Xp, the qq contribution is dominant with T,qq{xi)q{x2) ~ u{xi)u{x2) + d{xi)d{x2) 
and the change in the projectile sea quark distribution is less than that of the gluon 
distribution. When gg fusion dominates, a{xp) decreases. In the NRQGD model, 
the qg contribution tends to balance this difference, leading to the flatter a{xp) for 
Xp < 0, particularly at 800 GeV. As Xi;' approaches zero, the change in all the 
distributions becomes smaller. The change in xi due to the minimum BH loss is 
only ~ 20% at xp ~ 0, decreasing to less than 8% at xp = 0.1. Note also that 
the predicted minimum loss for gluons, Eq. (^), and the original BH loss, Eq. (^), 
result in a similar a{xp) at forward xp even though the A dependencies of the two 
are different. That is because the shift in xi is reduced at large xp, the change in all 
parton densities is small when Xi ~ 0.1 — 0.9 and the original BH model dE/dz is very 
similar to dEg/dz for the minimum BH loss. At 800 GeV, the drop at large is due 
to loss by valence quarks since at large xp, Egg(xi)g(x2) ~ u{x2)u{xi) + d{x2)d{xi). 
At 120 GeV, the effect of the loss is larger since Axi oc 1/s. The correspondingly 
higher x[ values at 120 GeV reduce the gluon distribution shift relative to qq. This, 
as well as the greater importance of qq annihilation, results in the different shapes of 
a(xp) for the two energies at negative xp. 

The Drell-Yan results are similar to the calculated ip results except that the Drell- 
Yan loss is weaker at larger xp. Part of the difference is because Xi and the scale 
M at which the parton densities are evaluated are both greater than for the ip. We 
can also see that the maximum BH loss is almost certainly too large to explain the 
current Drell-Yan results. Indeed, at 120 GeV, a{xp) barely appears on the plot. 
The Drell-Yan results are similar to the calculated ip results except 
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The large change in xi, appearing as large Axi, suggests that the calculation may 
not be apphcable for Axi > Xi. At 800 GeV, Axi = Xi occurs when 0.03 < Xi < 0.09, 
depending on the loss estimate, corresponding to a minimum xp of —0.1 < xp < 0.02. 
At 120 GeV, the xi values are larger, xi ~ 0.1 — 0.3 for the original and minimum BH 
loss estimates and 0.6 for the maximum BH estimate, corresponding to Xi? ~ -0.2-0 
and 0.5 respectively. In our calculations, we will apply the model over all xp- 

We finally note that neither of these initial state models of energy loss alone can 
reproduce the data. CM loss does not have the same curvature of the data at large 
xp while BH loss is too weak at large xp and too strong at low xp. 



7.2 Final State Loss 



The second model of energy loss we consider is applicable only to the quarkonium 
system and not to Drell-Yan production which does not involve color confinement in 
the final state 
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When a cc pair is produced in a color octet state, it has to emit 
a soft gluon in order to produce the final-state ip or This cc can propagate some 
distance, essentially longer than its path through the nucleus, before the soft gluon is 



finally emitted. This is because the Landau-Pomeranchuk-Migdal effect in QCD 
causes a delay in the emission of the third gluon to neutralize the color of the cc state 
due to successive interactions of the colored cc pair in the medium. However, each 
successive interaction of the cc pair degrades its momentum. 
This final-state loss model, 



developed by Kharzeev and Satz ||T^ and referred 
to as KS loss here, is applicable only when the cc pair interacts in the color octet 
state, essentially for xp > 0. After n interactions along its path length before leaving 
the target, the pair's momentum is reduced by ~ kLa where k, the hadronic string 
tension, is determined from lattice studies of confinement between colored objects, 
K ~ (9/4) GeV/fm [|l^, and La is the distance the pair has traveled through the 
target, calculated for the nuclear shape distributions in Ref. pO[. A ip state observed 



at a given xp has actually been produced with a higher value, xp/6, where 6 
1 — kLa/P-4, and is the ip momentum in the center of mass frame. 
The Xp distribution Ga{xp) then has two parts ||17|| , 



Ga{xp) oc SaG,{xp) + (1 - _ ^^/s) 





(59) 



where Gp{xp) is the xp distribution in pp interactions and Sa is the survival probabil- 
ity for the cc pair not to break up on its way out of the target, calculated in Eq. ( p9D 
for pure octet production. The second term includes the scatterings in the target 
that cause the shift in xp. The effect of Eq. (|^) does not produce an integrated 
"0 suppression: the integrated a in Eq. (|lD is unchanged with this mechanism, only 
a{xp) changes due to the shift in Xi;'. 

The resulting xp dependence is shown in Fig. O for Xi;' > 0. We have illustrated 



the effect for the MRST LO parton densities at 800 GeV with three different color 
octet cross sections: 1 mb, 20 mb and 40 mb. A 20 mb cross section was chosen 
originally to be as large as a typical meson-nucleon inelastic cross section at the 



same energy, a^hs ~ f^'jv ~ 20 mb. The 1 mb cross section shows a minimal effect 
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Figure 11: The A dependence of ip production assuming KS loss for xp > 0. Octet 
cross sections of 1 mb (solid), 20 mb (dashed) and 40 mb (dot-dashed) are calculated 
with the MRST LO parton densities at 800 GeV. At 120 GeV, a 40 mb octet cross 
section is assumed (dotted). 



while (Tabs ~ 40 mb sets the scale for the maximum effect since then (Tabs > c^pp • 
While a large octet cross section is needed to produce a strong effect at xp > 0.5, the 
normalization amplifies a{xF ^ 0) so that the shape of the dependence is significantly 
different from the behavior of the E866 data at low xp. Including shadowing will 
further increase a for xp ~ 0, as we will see. Due to its nature, this model is limited 
to the case where all ip^s are assumed to be produced in pure color octet states. 

It is clear that KS loss alone cannot account for the shape of the ip and ip' data 
in Fig. ^ A combination of effects is needed. 



8 Intrinsic Charm 

The wavefunction of a proton in QCD can be represented as a superposition of Fock 
state fluctuations, e.g. \uudg), \uudqq), \uudQQ), ...of the \uud) state. When the 
projectile scatters in the target, the coherence of the Fock components is broken and 
the fluctuations can hadronize |TB|, . These intrinsic QQ Fock states are dominated 
by configurations with equal rapidity constituents so that, unlike sea quarks generated 
from a single parton, the intrinsic heavy quarks carry a large fraction of the parent 
momentum |]18| . 



The frame-independent probability distribution of a 5-particle cc Fock state in 
the proton is 

- N^a.im)— — ^— rj^--^ , (60) 



dxi--- dx5 [ml - YH=i{fm 

where normalizes the \uudcc) probability, P;^. The delta function conserves lon- 
gitudinal momentum. The denominator is minimized when the heaviest constituents 
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carry the largest fraction of the longitudinal momentum, (xg) > (xg), maximizing 
Pt We choose m„ = 0.45 GeV and = 1-8 GeV 12011. 



The intrinsic charm production cross section from the 5-particle state can be 
related to P;^ and the inelastic pN cross section by 

<{pN) = Pl^'^M^^ ■ (61) 

The factor of /i^/4m^ arises from the soft interaction which breaks the coherence of 
the Fock state. We assume that the NA3 diffractive ip cross section [^], the second 
term in Eq. (0) proportional to A^^ , can be attributed to intrinsic charm and find 
/i^ ~ 0.1 GeV2. 

While the total intrinsic charm cross section is relatively easy to define, there 
are some uncertainties in the relative weights of open charm and ip production from 
an intrinsic charm state. In general, the tp production cross section is significantly 
smaller than the open charm production cross section. There are several factors that 
can suppress ip production relative to open charm in standard charmonium production 
models such as the GEM and NRQGD as well as in the intrinsic charm model. As 
in the GEM, the probability to produce a ip from an intrinsic cc state is proportional 
to the fraction of intrinsic cc production below the DD threshold. The fraction of cc 
pairs with 2mc < m < 2m/) is 

Mml, dp. I rs dp. 

^-/^ = L / L ^ ' (^2) 



typically smaller than that obtained in the GEM |^T|. However, as discussed in 
Section 2.1, not all cc pairs below the DD threshold will produce a final state tp. The 
fraction that actually become ip^s is rather small, on the order of 2.5% in the GEM 



2T[| . Since the additional suppression factors involved in the intrinsic charm model 



are not completely fixed rather than discuss all the uncertainties here, we will 



use an effective intrinsic charm probability, P^^ ■ The EMG charm structure function 
data is consistent with Pj^ = 0.31% for low energy virtual photons but P;^ could be 
as large as 1% for the highest virtual photon energies p6| , |67[ . Typically the more 
conservative result is used but in this paper, we will use the larger value in most 
of our calculations and show the effect of reducing and/or eliminating the intrinsic 
charm component. 

Including a delta function to combine the Xc and Xc in the ip state, the ip xp 
distribution from intrinsic charm is 

S = / n d^^, ^(^^ - - • (63) 

Only the 5 particle Fock state is considered. The intrinsic charm contribution is 
included as in Eq. (|^) with (3 = 0.71. The total A dependence for intrinsic charm 
alone is shown in Fig. |l^ assuming P^^ = 1% and 0.31% with both charmonium 
production models. The contribution is symmetric around xp = since the projectile 
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Figure 12: The A dependence of intrinsic charm at 800 GeV (sohd) and 120 GeV 
(dashed). In (a) and (c) an effective production probabihty of 1% is assumed in the 
GEM and in NRQGD respectively while in (c) and (d) i^f = 0.31% is assumed in 
the GEM and in NRQGD. 
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and target fragmentation regions are treated equally. Figure [T^(a) and (c) show a{xp) 
in the CEM. A larger effect is seen at high xp with the NRQCD model, Fig. 0(b), 
because the NRQCD xp distribution is narrower. Since charmonium production 
models outlined in Section 2 have a larger energy dependence than the intrinsic charm 
cross section in Eq. (pT|), intrinsic charm is more important at 120 GeV. When P-^^ = 
0.31%, the intrinsic charm contribution to the total A dependence is quite small and 
only significant for the largest xp values due to the reduced A dependence of the 
mechanism. Assuming a 1% probability enhances the intrinsic charm effect at large 
Xp and even suggests that intrinsic charm can influence a{xp) at xp ~ 0. 



9 Results and Predictions 

We now have a comprehensive model with which we can confront the nuclear de- 
pendence of ipy 4''^ cind Drell-Yan production. The nuclear effects included in the 
model are shadowing of the parton distributions, energy loss, nuclear absorption, co- 
mover interactions, and the 'diffractive' intrinsic charm component. It is clear from 
an examination of the individual nuclear effects shown in Figs. that no single 
mechanism can correctly predict the shape of the E866 xp data. 

To compare with the preliminary data, we calculate a{xp) of ip and ip' production 
with all three models of energy loss and all shadowing parameterizations. For pure 
octet absorption, we use all three models of energy loss. Only GM and BH loss are 
used with pure singlet absorption and the combination of singlet and octet produc- 
tion. The CEM model with the MRST LO parton distributions is used to calculate 
charmonium production for pure octet and pure singlet absorption. NRQCD is used 
as the basic production model for the octet/singlet combination, Eq. (pSj). All three 
shadowing parameterizations are used in each case. We use an effective intrinsic 
charm probability of 1% but will examine the relative importance of intrinsic charm 
to the overall description of the large xp E866 data. The absorption cross sections 
are chosen so that the shadowing parameterization gives reasonable agreement with 
the magnitude of a{xp) for both ip and ijj' production ai xp > with GM loss. We do 
not actually make detailed fits to the data to obtain the cross sections. The resulting 
absorption cross sections are given in Table |l| for GM and BH loss. The KS model 
of energy loss is always calculated with an octet absorption cross section of 40 mb. 
The NA3 Pt/^H ratio as a function of xp at 200 GeV is also compared to the model 
calculations. We make predictions of the ip and ip' A dependence at 120 GeV. Finally, 
we show the combined effects of shadowing and initial-state energy loss on Drell-Yan 
production at 800 and 120 GeV. 

We first compare our full model results with the preliminary E866 ip data in 



Figs. |T3|, |1J, and |T5|. Each figure shows the difference in the shadowing mechanisms 
for each type of energy loss with a particular absorption mechanism. In Fig. |T3|, the 
pure octet absorption mechanism is shown. The results with GM loss. Fig. |l3|(a), 
best reproduce the general trend of the data for xp > 0.1. In general the agreement 
is worse a.t xp < because X2 is in the antishadowing region where antishadowing 
of gluons enhances a{xp), see Figs. Bl-Rl The KS loss model is typically above the 
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Figure 13: All effects are compared with the preliminary E866 data assuming 
pure octet absorption. In (a) and (b), GM and KS loss are assumed. Energy loss 
effects associated with the BH bound are shown in (c), (d), and (e) for the estimated 
maximum and minimum loss and the original bound respectively. All calculations are 
in the GEM with the MRST LO parton densities. The curves represent shadowing 
with the 5*3 (solid), 5*2 (dashed) and (dot-dashed) parameterizations. 
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data, except for the 5*2 parameterization, and is inconsistent with the shape of the 
prehminary E866 data at xp < 0.2. Since the KS model is only applicable for xp > 0, 
we make no further calculations with this model. The calculations of a{xF) with BH 
loss. Fig. |l3](c)-(e), do not match the data very well at low xp, particularly for the 
maximum BH loss. While the results with the maximum loss, Eq. (p6D, produce the 
largest reduction at large xp, the negative xp region is far off due to the drop in a 
at negative xp. The curvature of a{xp) changes at xp ~ 0.1 — 0.25, the point at 
which the slope of the BH loss flattens in Fig. ^OKa). Better results are achieved with 
the lower estimates of the BH loss, the minimum estimate, Eqs. ( p7D and and 
the original suggestion, Eq. (0). The data are somewhat overestimated at xp < 0.2 
with the minimum loss but the overestimate is slight for the original BH loss with 
the A^^^ dependence. Choosing a smaller absorption cross section would improve 
the agreement with the data at low xp although it would worsen the agreement at 
Xp > 0.25. We also note that none of the absorption cross sections are greater than 3 
mb, already more than a factor of two less than the 7.3 mb effective absorption cross 
section found in Ref. 



Pure singlet absorption, shown in Fig. |14|, results in somewhat poorer agreement 
with the data than pure octet absorption because the ip is always produced outside 
the target when > at 800 GeV. Therefore changing the absorption cross section 
would not improve the agreement with the data. The choice of parton distribution 
function in the CEM model results in small changes in the shape of a{xp) and does 
not influence the overall agreement with the data. 

A combination of octet and singlet absorption in the NRQCD ip production model 
produces rather good agreement for all shadowing parameterizations when the GM 
model is applied. Fig. |15|(a). The difference in curvature a.t xp > between the 
calculations with BH loss and the preliminary data are larger than in the CEM 
because the A dependence of shadowing in NRQCD is weaker at positive xp than 
that of the CEM, see Figs. ^[|^. As explained in Section 6, the difference in the A 
dependence of the two production models is due to the chosen x values and the charm 
quark mass which sets the scale for evolution. 

In all cases, the most striking disagreement of the GM loss model with the data 
occurs at xp < 0.1 when the calculated a slightly overshoots the data due to the anti- 
shadowing of the gluon distribution. The shape of a{xp) here depends most strongly 
on the shadowing parameterization since the other xp dependent contributions are 
rather slowly varying. None of these parameterizations produce the same curvature 
as the data and, even if they did, the additional absorption required for the calcula- 
tions to agree with the data would ruin the agreement of the model with the data at 
forward xp. Increasing absorption a.t xp < hj artificially enhancing the co mover 
density would not significantly improve the agreement. 

The pure octet and pure singlet calculations with BH loss are in reasonably good 
agreement with the data for xp > 0.2 The maximum estimated loss is in clear dis- 
agreement with the data at allxp, both in shape and in magnitude. This disagreement 
would persist, even if final-state absorption of the ip were ignoredPl. It is also then 



^S. Gavin has addressed the E772 -0 data with a combination of BH-type loss and shadowing 
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Figure 14: All effects are compared with the preliminary E866 ip data assuming 
pure singlet absorption. In (a), GM loss is assumed. Energy loss effects associated 
with the BH bound are shown in (b), (c), and (d) for the estimated maximum and 
minimum loss and the original bound respectively. All calculations are in the GEM 
with the MRST LO parton densities. The curves represent shadowing with the 5*3 
(solid), 5*2 (dashed) and Si (dot-dashed) parameterizations. 
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Figure 15: All effects are compared with the preliminary E866 ip data assuming 
a combination of octet and singlet production and absorption. In (a), GM loss is 
assumed. Energy loss effects associated with the BH bound are shown in (b), (c), 
and (d) for the estimated maximum and minimum loss and the original bound re- 
spectively. All calculations are in NRQCD with the CTEQ 3L parton densities. The 
curves represent shadowing with the 5*3 (solid), 5*2 (dashed) and 5*1 (dot-dashed) 
par ameterizations . 
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Figure 16: The effective probability of intrinsic charm is varied for pure octet produc- 
tion with (a) GM loss and (b) the minimum BH loss. The curves represent an effective 
intrinsic charm probability of 1% (sohd), 0.31% (dashed) and 0% (dot-dashed). 

unlikely to produce results consistent with the minimal A dependence of Drell-Yan 
production, especially at lower energies, as we discuss later. The agreement with 
the minimum and original BH calculations is reasonable for xp > 0.2. At lower xp 
values, the change in Axi due to the energy loss is large. However, since Axi ~ Xi 
at Xp ~ 0, the model is at the limit of applicability and therefore the magnitude of 
the disagreement is suspect. 

To show the influence of intrinsic charm, we take the 5*3 shadowing parameteriza- 
tion with GM and minimum BH loss and vary P-^^ between and 1% in Fig. We 
choose pure octet production and absorption because the agreement with the data 
seems to be among the best, see Figs. . Since the GM loss mechanism alone 

causes strong reduction in a at large xp, see Fig. |^, including intrinsic charm does 
not have a large effect. It would appear from Fig. |l6|(a) that P^^ = 0.31% agrees 
best with the data although the agreement is reasonable in all three cases. The same 
is true for the combination model but pure singlet absorption would require a larger 
intrinsic charm probability to agree with the data. On the other hand, the relatively 
good agreement of the minimum BH loss calculations with the data at large xp is 
due to the intrinsic charm contribution. Without intrinsic charm with P^^ = 1%, the 
model calculations would not agree with the data. The minimum and original BH 
loss models affect a weakly at positive xp and shadowing alone can only reduce a 
to ~ 0.85 as Xp ^ 1 with the 5*2 parameterization, see Fig. ^ Thus increasing the 
relative intrinsic charm contribution is the only way to produce agreement with the 
data at large xp. This is clearly shown in Fig. |16|(b). With no intrinsic charm, a{xp) 
is relatively fiat at large xp. Similar results are obtained with the pure singlet and 
combination absorption models. Note that for both loss mechanisms, intrinsic charm 

without final-state effects and found a similar level of agreement as is seen here . 
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only affects the shape of a{xp) at xp > 0.25. 

The corresponding ip' calculations are compared to the data in Fig. |l^ with the 
5*3 shadowing parameterization. For each absorption model, we show the A depen- 
dence with the GM loss and the three estimates of BH loss. Because of the greater 
uncertainties in the data, none of the calculations are fully incompatible with the 
data at forward xp. The largest discrepancies between model and data are at low 
to negative xp. GM loss produces the largest a because of negligible loss effects at 
Xp ^ combined with antishadowing. The minimum and original BH loss models 
agree relatively well with all the data for the three absorption models. Note that 
these two loss models coincide at large xp both because of their similar behavior at 
large xp, see Fig. and the intrinsic charm contribution at forward xp. All the 
calculations result in a slightly lower a for the ip' than the ip due to the larger comover 
cross section, a^'co- 

The A dependence of and tp' production has been shown to be similar in previous 
measurements albeit not to high precision. To compare the two results here, we 
calculate the integrated a in the interval —0.2 < xp < 0.8 for all energy loss models 
and all shadowing parameterizations. The results are shown in Tables The 
change in a between ip and ip' at 800 GeV is small, typically a 2 — 3% difference. One 
might expect that for pure octet absorption, the integrated a would be identical for 
ijj and Ip'. However, the pA comover interactions are treated assuming formed ip and 
ijj' interact with secondaries and the ip' comover absorption cross section is larger. 
Thus, even though the comover interaction cross sections are typically significantly 
smaller than the corresponding nucleon absorption cross sections, the difference is 
large enough to cause the observed 2% shift in the integrated a in the Tables. (See 
also Fig. ^(d) which highlights the differences in the assumed comover cross sections.) 
Indeed, without this difference in the cross sections, the model calculations would not 
agree as well with the data in Fig. [T^ 

We note that we do not expect our values of a to agree in detail with the integrated 
data because our estimates do not agree with the preliminary data at all xp. The 
GM model always overestimates the data at low to negative xp. Thus the GM results 
can then be expected to overestimate the integrated a of the data. Typical a values 
are between 0.94 and 0.98. The original and minimum estimates of BH loss typically 
underestimates the low and negative xp data and should therefore underestimate the 
measured total a. In this case, 0.87 < a < 0.94, similar to a = 0.91 |Q. Even though 
the nuclear absorption cross sections are small, the effective absorption can be large, 
compatible with that obtained assuming absorption is the only source of the ip A 
dependence. 

The integrated a also depends on the shadowing parameterization. Typically 
a with 5*2 is largest with the octet absorption because the 5*2 parameterization is 
not available for A = 9 and thus treats the Be nucleus like a proton. Therefore, 
although the 5*3 parameterization has a larger gluon antishadowing effect than 5*2, the 
calculated value of a is larger for 5*2 near xp ^ 0. This difference, along with smaller 
absorption cross sections used with 5*2, see Table |T], results in a larger integrated a 
with S2 for pure octet absorption. Table 0. Since the ip and ip' are formed after they 
have left the nucleus in the pure singlet case at 800 GeV, the a obtained with S3 
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Figure 17: All effects are combined and compared with the preliminary E866 ip' data 
1^. In (a), pure octet absorption is assumed. The results for pure singlet absorption 
are given in (b) and combined octet/singlet absorption in (c). All calculations in (a) 
and (b) are in the CEM with the MRST LO parton densities while the calculations 
in (c) are in NRQCD with the CTEQ 3L parton densities. The parameterization 
is used for calculations with different energy loss models. The GM loss is shown in 
the solid curve while the dashed, dot-dashed, and dotted curves are calculations with 
the maximum, minimum, and original BH loss estimates respectively. 
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is larger at this energy, see Table 0. In the combination octet/singlet calculations, 
Table ^, a is very similar for 5*2 and 5*3 . 

A measurement of the A dependence at 120 GeV has been proposed. Such 
an energy would be the closest to the NA50 Pb+Pb measurement. The most precise 
lower energy data with a proton beam was taken at 200 GeV by the NA3 collaboration 
0. We compare our calculations of the ratio Pt/^H to their data in Fig. 18. The 
three absorption mechanisms are shown with the data for GM loss and the three BH 
loss estimates with the 6*3 shadowing parameterization. We see that essentially none 
of the calculations contradict the large xp data, presumably due to the relatively poor 
statistics for xp> 0.5. The best agreement at all xp vs, obtained with minimum and 
original BH loss with the pure singlet and combination absorption models. GM loss 
tends to underestimate the low xp data. The maximum BH loss estimate produces 
a higher ratio at low xp^ as seen in Fig. |l8|(a) and (c), far above the data. On the 
whole, the results at 200 GeV concur with those at 800 GeV. 

Figures |l^ and ^ show the predictions for the '0 and ^' A dependence at 120 
GeV. The calculations with GM loss show a plateau-like behavior at this energy 
and are very similar to each other over all xp. This is due to the widening of the 
gluon antishadowing region over the xp interval, see Figs. In the pure singlet 
model, since the i\} and ■?/'' may be produced inside the target at forward xp^ there 
is a peak in the calculated il) a{xp^ at xp ~ 0.2, see Fig. |l9|(b). This peak is 
shifted slightly forward for production since the ■?/'' singlet absorption cross section 
is larger, see Fig. pO|(b). There is no forward xp peak when the octet model is 
considered because the nucleon absorption is treated identically at all xp. Since the 
combination octet /singlet model includes both types of absorption, there is a rather 
wide plateau over —0.2 < xp < 0.3. The shape of a{xp^ at 120 GeV could therefore 
help distinguish between absorption models. However, BH loss yields very different 
expectations at the lower energy. At low xp^ the BH loss calculations are governed 



by the decrease shown in Fig. |T0|. This effect is enhanced with singlet absorption 
since the ■?/' and ■?/'' have some probability to interact with their full cross sections. 
The maximum BH loss is so large that the first component of Eq. (^) is less than 
the intrinsic charm contribution so that rather than having a peak in a at low Xi?, a 
minimum is seen instead. Results with the GM model and the minimum and original 
BH losses are similar for xp > 0.5 due to intrinsic charm. 

Note that at this energy, the integrated values of a, also given in Tables 0-^, 
indicate that the differences between the ijj and ijj' A dependence are again on the 2% 
level except when pure singlet absorption is considered. In the pure singlet case, the 
absorption cross section grows more slowly with at 120 GeV and some portion of 
the ijj and ip' resonances are formed inside the target a.t xp > 0. Even though the ip' 
formation time is larger than the ip formation time, the singlet absorption is larger 
for the ip' due to its increased cross section. High statistics measurements of the ip 
and ip' A dependence at this energy would set bounds on the importance of singlet 
production of the resonances since there is a 3% difference in a between ip and ip' 
at 120 GeV in the singlet case. Note also that the ip and ijj' octet and combination 
A dependencies are within 1-2% of each other at this energy, generally a smaller 
difference than at 800 GeV, because the comover density is lower. 



42 



Pu 
b 

\ 

M 

b 





0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 



-P 

b 

\ 

Oi 

b 
13 




Xx 



1.00 



b 

\ 
b 



Figure 18: All effects are combined and compared with the NA3 ip data 0]. In (a), 
pure octet absorption is assumed. The results for pure singlet absorption are given 
in (b) and combined octet/singlet absorption in (c). All calculations in (a) and (b) 
are in the CEM with the MRST LO parton densities while the calculations in (c) 
are in NRQCD with the CTEQ 3L parton densities. The 5*3 parameterization is 
used for calculations with different energy loss models. The GM loss is shown in the 
solid curve while the dashed, dot-dashed, and dotted curves are calculations with the 
maximum, minimum, and original BH loss estimates respectively. 
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Figure 19: Predictions are made for ip production at 120 GeV. In (a), pure octet 
absorption is assumed. The results for pure singlet absorption are given in (b) and 
combined octet/singlet absorption in (c). All calculations in (a) and (b) are in the 
GEM with the MRST LO parton densities while the calculations in (c) are in NRQCD 
with the CTEQ 3L parton densities. The parameterization is used for calculations 
with different energy loss models. The GM loss is shown in the solid curve while the 
dashed, dot-dashed, and dotted curves are calculations with the maximum, minimum, 
and original BH loss estimates respectively. 
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Figure 20: Predictions are made for ip' A dependence at 120 GeV. In (a), pure octet 
absorption is assumed. The results for pure singlet absorption are given in (b) and 
combined octet/singlet absorption in (c). All calculations in (a) and (b) are in the 
GEM with the MRST LO parton densities while the calculations in (c) are in NRQCD 
with the CTEQ 3L parton densities. The parameterization is used for calculations 
with different energy loss models. The GM loss is shown in the solid curve while the 
dashed, dot-dashed, and dotted curves are calculations with the maximum, minimum, 
and original BH loss estimates respectively. 
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It is typically assumed that the energy dependence of a is small. This is true for 
both the GM loss and the minimum and original BH loss estimates where a changes 
between 2 and 6%. However, the energy dependence of the maximum BH loss is much 
stronger. There is a 10-25% decrease of the integrated a with this estimate for the 
pure octet and singlet absorption mechanisms, as can be seen in Tables ^ and ^ The 
energy dependence of the combination absorption model is typically stronger for all 
BH loss estimates because there is a stronger energy dependence of BH loss when ip 
production is calculated in NRQCD, compare Figs, ^(a) and (b). As seen in Table ||, 
the resulting energy difference in a is as large as 35% for the maximum BH estimate 
but is only 3-7% for the lower BH loss estimates. Precision measurements of a at 
800 and 120 GeV would help eliminate models. For example, the maximum BH loss 
could be ruled out but discerning the difference between the A^^^ dependence of the 
loss in the original BH estimate and the A'^^^ dependence of the minimum BH loss, 
Eq. (p7D, could be difficult. 

Finally we compare the difference between the Drell-Yan A dependencies with 
shadowing and energy loss at 800 and 120 GeV in Figs. ^ and ^ respectively. Recall 
that the Drell-Yan mass range is 4 < M < 9 GeV. All three shadowing param- 
eterizations are shown for each energy loss estimate. We have included the E772 
measurement of the Drell-Yan A dependence based on D, Ca, Fe and W targets p[ 



in Fig. |2T]. Since the GM loss parameter, eg, is tuned to this data without shad- 
owing, the full model calculation overestimates the A dependence at large xp- All 
the calculations with BH loss predict a significantly less than linear A dependence 



for Xi? < at 800 GeV, as shown in Fig. 10. The maximum BH loss estimate is in 



complete disagreement with the E772 data which shows a rather mild A dependence. 
The deviation of the model from the data in Fig. pT|(b) suggests that the maximum 
estimated BH loss is certainly too large. Better agreement with the data is found 
with the original BH loss in Fig. ^(d) for the 5*1 and 5*3 shadowing parameteriza- 
tions although the calculations consistently overestimate the A dependence. When 
the difference in the color factors is included, the minimum BH quark loss, Eq. (^), 
agrees reasonably well with the data. Indeed, this model is the only one that follows 
the trend of the E772 data over all xp- 

The 120 GeV calculation shows a significant energy dependence of a when BH 
loss is considered. The A dependence is now calculated for W and Be targets, as 
used in the E866 experiment. A precision measurement of at this energy could 

reveal if the A dependence is less than linear here. There is, up to now, no data 
on the A dependence of Drell-Yan production with a proton projectile below 200 
GeV. A measurement of a at 200 GeV is for low masses, 1.7 < M < 2.7 GeV ||69 |. 



Earlier measurements at lower energies are only available with pion beams [^, |7T 
The lowest energy pion beam data, at 140 GeV, obtained a = 0.980 ± 0.006 ± 
0.013 which leaves room for a less than linear A dependence at the lower proton 
beam energy. Other pion data below 280 GeV does not have enough statistical 
significance to determine whether a deviates from unity at this energy, particularly 
at large xp- 

Since the E866 Drell-Yan data seems to indicate that only shadowing is neces- 
sary to explain the Drell-Yan A dependence, in Fig. Pl(a), we compare the E772 A 
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Figure 21: Shadowing and energy loss in Drell-Yan production are combined and 
the A dependence calculated for W and D targets, shown at 800 GeV. The E772 
Drell-Yan A dependence is also shown. In (a), GM loss is assumed. Energy loss 
effects associated with the BH bound are shown in (b), (c), and (d) for the estimated 
maximum and minimum loss and the original bound respectively. All calculations 
are with the MRST LO parton densities. The curves represent shadowing with the 
S3 (solid), 5*2 (dashed) and 5*1 (dot-dashed) parameterizations. 
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Figure 22: Shadowing and energy loss in Drell-Yan production are combined and 
predictions of the A dependence for W and Be targets given at 120 GeV. In (a), 
GM loss is assumed. Energy loss effects associated with the BH bound are shown 
in (b), (c), and (d) for the estimated maximum and minimum loss and the original 
bound respectively. All calculations are with the MRST LO parton densities. The 
curves represent shadowing with the S3 (solid), S2 (dashed) and 5*1 (dot-dashed) 
parameterizations. 
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dependence with calculations of shadowing effects alone. All three shadowing pa- 
rameterizations are in reasonable agreement with the data. That should be expected 
because the 5*2 and 5*3 parameterizations included these data in their fits. A com- 
parison of Figs. |2l|(c) and ^(a) shows that the minimum BH quark loss has a very 
weak effect on the shape of a{xF). The energy loss is only obvious in the curves at 
Xp < 0.1 and xp > 0.85 and has little apparent influence on the agreement with the 
data. Thus, at least for this case, energy loss may be present but the effects nearly 
indistinguishable from those due to shadowing alone. The 120 GeV predictions with 
shadowing only are shown in Fig. ^(b). Contrary to the results in Fig. shadowing 
alone predicts a negligible influence on the A dependence at the lower energy. A clear 
distinction can be made between models which include energy loss and those which 
do not with this measurement. Even a very small quark loss, such as that in the min- 
imum BH quark loss estimate, causes the A dependence to be less than linear at 120 
GeV while shadowing alone would suggest that the A dependence is either exactly 
linear or slightly greater than linear over all xp- A high statistics measurement of the 
Drell-Yan A dependence at 120 GeV could decisively settle the issue. 
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Figure 23: Shadowing effects on Drell-Yan production at 800 (a) and 120 (b) GeV. All 
calculations are with the MRST LO parton densities. The curves represent shadowing 
with the 5*3 (solid), 5*2 (dashed) and 5*1 (dot-dashed) parameterizations. Note that 
in (a), a is calculated with W and D targets for comparison with the E772 data 
while the calculations at 120 GeV assume W and Be targets. 



The integrated a values for all the calculations at 800 GeV and 120 GeV assuming 
W and Be targets as in the E866 experiment are given in Table |^. The choice of Be 
as the lowest mass target results in lower values of a than for a very light target like 
D since nuclear effects are larger in Be than D. When the GM loss is considered, all 
shadowing parameterizations are consistent with a 1% determination of a = 1 at 800 
GeV. At the lower energy, a is up to 2% away from unity with the 5*3 parameterization. 
These results can be expected from the shadowing of the sea quark distributions in 
the 'antishadowing' range, 0.1 < a; < 0.3, in the 5*3 parameterization. (Compare 
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Figs, ^(c) and |^(c) with Fig. ^(d).) If the original BH model is correct, a less than 
linear integrated A dependence should already have been observed at 800 GeV. Since 
this is not the case, it seems that all but the minimum estimates are too large to 
explain the current Drell-Yan data. We note that even a reduced loss which agreed 
with the 800 GeV data would still predict a < 1 at 120 GeV due to the energy 
dependence of the BH model which should be observable with high precision data. 
These results can be compared to the a predicted for shadowing alone, all compatible 
with unity at 120 GeV. 



10 Conclusions 

In this paper we have tested a number of different nuclear effects: nuclear absorp- 
tion, comover scattering, nuclear shadowing, energy loss and the A dependence of 
intrinsic charm. We have not attempted to make precise fits but rather check the 
resulting shape of a{xp) for a number of processes that are expected to contribute 
to ip 'suppression' in pA collisions. The processes with the strongest influence on the 
xp dependence are shadowing and energy loss. 

It is clear from the comparisons of the ip and ip' calculations with the data that 
a single mechanism cannot describe the shape of a{xF) for all xp- Combining all 
effects can explain some portion of the data, depending on which model of energy 
loss is assumed. A constant energy loss, a la Gavin and Milana [jl3l, can describe 



the forward data when combined with the other effects discussed here but results in 
values of a too large at low xp due to gluon antishadowing. 

An energy dependent model of energy loss, like that proposed by Brodsky and 
Hoyer with refinements by Baier et al. |5^, |6^, is less influenced by antishadowing 



because the dependence Axi oc 1/xiS produces a strong xi shift at negative xp, even 
strong enough to counteract the gluon antishadowing. At xp < 0, the application 
of the model becomes problematic because Axi grows larger than Xi. Therefore the 
results in this region should be treated cautiously. This type of model alone fails to 
explain the data at larger xp because the xi shift becomes too small to cause a large 
enough change in the parton distributions for xp > 0.25. The BH loss estimates are 
rather crude. For example, the minimum BH loss was calculated assuming that the 
gluon distribution in the proton remains relatively constant over the Xi range of the 
data. If the nuclear gluon density as a function of xi was considered instead, the 
results may be more compatible with the data. It is unlikely however that such an 
effect could be significant at large xp since the gluon distribution decreases strongly 
at large Xi, thereby weakening the effect at large xp. Thus if the BH model of energy 
loss is correct, a combination of BH loss and shadowing alone cannot describe the 
data, further strong absorption at large xp is still needed. 

The apparently stronger absorption at large xp seen in the NA3 data 0] was the 
motivation for introducing ip production by hadronization of intrinsic charm states 
20]. Indeed, without intrinsic charm, a{xp) at large xp would be even further above 



the data, see Fig. 16. We have used an effective intrinsic charm production probability 



of 1%, within the uncertainties of the production probability determined from a fit 
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to deep-inelastic scattering . 

We have primarily used the MRST LO and CTEQ 3L parton densities in our 
calculations. While some deviations in the expectations of a{xF) appear for other sets 
of parton densities, these are not significant enough to change our general conclusions. 
The most important effect in the determination of the A dependence is the nature of 
the energy loss. Understanding this loss requires correlation of the ip A dependence 
with that of other processes. 

The preliminary data thus seem to suggest that final-state absorption, regardless of 
the mechanism, is not as strong as previously expected from studies of the integrated 
A dependence |Tl|] when other nuclear effects are included. A smaller cross section is 
needed than determined from absorption alone. This would be true even if comover 
effects are neglected since they are very small in pA interactions. If the energy loss of 
quarks and gluons is treated on an equal footing, the BH loss mechanism results in a 
stronger A dependence than required for the Drell-Yan data. Indeed, the minimum 
BH loss shows that treating the energy loss of quarks and gluons separately can lead 
to qualitative agreement with the ip and Drell-Yan A dependence. 

Further data on the A dependence at 120 GeV could clarify the relative importance 
of octet and singlet states in the production and absorption of the charmonium states. 
A Drell-Yan measurement at this energy may decisively determine the importance 
of energy loss by the initial partons. In addition, precision measurements of the 
xp dependent absolute cross sections in pp collisions could show whether the xp 
distribution is closer to that expected from the CEM or NRQCD. 

The lower energy data can also provide an additional point of comparison to the 
NA50 measurements of ip suppression in heavy-ion collisions at the CERN SPS 0. 
Only nuclear absorption and some comover scattering has been used to compare to 
the XiT-integrated data |Tl|. Shadowing should also be included in the analysis 

7^. Interestingly, the a^xp) extracted from the E866 800 GeV data in the NA50 



xp region is larger than that obtained by NA50 between 158 GeV and 450 GeV 
An independent measurement at a similar energy could be very valuable. 
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octet 


singlet 


combination 




(^In (mb) 


cT^jv (mb) 


(T^T' (mb) 


</;']v'''* (mb) 


Si 


2 


8 


3 


2 


S2 


1 


5 


1 


1 


Ss 


3 


10 


3 


5 



Table 1: The ijj absorption cross sections used with each shadowing parameterization. 
Note that the corresponding ip' absorption cross sections are the same as those for 
the ip in the octet case and a factor of 3.7 larger for singlet production. In all cases 
the comover cross sections are a^co — 0.67 mb and cr^'co = 3.7(T^co- 



55 









^' 


Model 


S 


800 GcV 


120 GeV 


800 GeV 


120 GeV 






0.938 


0.922 


0.918 


0.907 


GM 




0.959 


0.953 


0.939 


0.938 






0.951 


0.922 


0.931 


0.907 






0.984 


0.985 


0.982 


0.979 


KS 


So 


0.980 


0.998 


0.978 


0.992 




& 


1.014 


1.002 


1.012 


0.996 






805 


661 


786 


648 




So 


0.823 


0.685 


0.803 


0.672 




S3 


0.817 


0.658 


0.798 


0.646 




Si 


0.872 


0.843 


0.852 


0.828 


min BH 


S2 


0.891 


0.872 


0.871 


0.857 




S3 


0.884 


0.842 


0.865 


0.827 




Si 


0.890 


0.850 


0.870 


0.836 


orig BH 


S2 


0.908 


0.875 


0.888 


0.864 




S3 


0.903 


0.851 


0.883 


0.836 



Table 2: The integrated value of a for ip and ip' production at 800 GeV and 120 GeV 
assuming pure octet absorption. Note that a is integrated over the range —0.2 < 
xp < 0.8 at both energies. 











Model 


S 


800 GeV 


120 GeV 


800 GeV 


120 GeV 




Si 


0.958 


0.916 


0.933 


0.882 


GM 


S2 


0.968 


0.944 


0.944 


0.914 




S3 


0.983 


0.924 


0.957 


0.886 




Si 


0.829 


0.658 


0.805 


0.631 


max BH 


S2 


0.833 


0.679 


0.811 


0.655 




S3 


0.854 


0.661 


0.830 


0.632 




Si 


0.894 


0.840 


0.870 


0.808 


min BH 


S2 


0.901 


0.865 


0.878 


0.836 




S3 


0.920 


0.846 


0.895 


0.811 




Si 


0.912 


0.848 


0.888 


0.816 


orig BH 


S2 


0.918 


0.873 


0.895 


0.845 




5';-! 


0.938 


0.855 


0.913 


0.821 



Table 3: The integrated value of a for ip and tp' production at 800 GeV and 120 
GeV assuming pure singlet absorption. Note that a is integrated over the range 
—0.2 < xp < 0.8 at both energies. 
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i\io(lc>i 


c 
o 


800 GoV 


120 Gov 


800 GeV 


120 Gov 




o 


0.948 


0.939 


0.926 


0.917 


GM 


'52 


u.y ^ ^ 


n no o 
0.983 


0.9OO 


0.9d3 




c 


0.978 


0.956 


0.954 


0.929 






0.804 


0.597 


0.788 


0.583 


T> TT 

max BH 


5^2 


U.c5oi 


U.Do / 


U.oi4 


U.DZO 




Q 

•J 3 


0.835 


0.612 


0.816 


0.593 




5l 


0.880 


0.830 


0.863 


0.817 


mill BH 


^2 


0.907 


0.872 


0.890 


0.860 




^3 


0.911 


0.847 


0.892 


0.830 




Si 


0.899 


0.853 


0.881 


0.835 


orig BH 


S2 


0.926 


0.896 


0.908 


0.879 






0.931 


0.871 


0.910 


0.847 



Table 4: The integrated value of a for and production at 800 GeV and 120 GeV 
assuming a combination of octet and singlet absorption. Note that a is integrated 
over the range —0.2 <xp< 0.8 at both energies. 



Model 


S 


800 GeV 


120 GeV 




Si 


0.998 


0.994 


GM 


S2 


0.997 


0.990 




S3 


0.993 


0.980 




Si 


0.916 


0.757 


max BH 


S2 


0.914 


0.753 




S3 


0.910 


0.740 




Si 


0.980 


0.967 


min BH 


S2 


0.979 


0.964 




S3 


0.974 


0.953 




Si 


0.973 


0.944 


orig BH 


S2 


0.972 


0.941 




S3 


0.968 


0.930 




Si 


0.997 


1.002 


no loss 


S2 


0.997 


0.999 




S3 


0.993 


0.989 



Table 5: The integrated value of a for Drell-Yan production at 800 GeV and 120 GeV 
for different model assumptions. Note that a is integrated over the range —0.2 < xp < 
0.8 at both energies. 
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